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Abstract. 

Calculations of the monoenergetic radial diffusion coefficients are presented 
for several configurations of the TJ-ff stellarator usually explored in operation. 
The neoclassical radial fluxes and the ambipolar electric field for the standard 
configuration are then studied for three different collisionality regimes, obtaining 
precise results in all cases. 



1. Introduction 

Radial electric fields are recognized to play a key role in the radial transport of stellara- 
tors. From the neoclassical transport point of view, they affect the particle orbits [1, 2]: 
for low-collisionality plasmas, they suppress the unfavorable \ jv regime [3] and allow 
for the formation of electron transport barriers, see Ref. [4] and references therein. 
Additionally, radial electric fields and plasma rotation are tightly connected: it is 
considered that sheared E x B flows are likely to reduce the edge turbulence level 
thus facilitating access to High confinement (H) mode, see Refs. [5, 6] and references 
therein. A review on internal transport barriers and H mode in helical systems can be 
found in Ref. [7]. Both effects have been measured at the heliac TJ-II [8]: transitions 
to core electron root confinement have been observed in Electron Cyclotron Heated 
(ECH) plasmas [9]. Transitions to H mode have been documented [10, 11] together 
with mean and low frequency oscillating sheared E x B flows in plasmas heated by 
Neutral Beam Injection (NBI). 

Neoclassical transport theory allows to predict the radial electric field in helical 
devices by means of the ambipolarity condition. Examples of these calculations exist 
for devices such as W7-AS [12] and LHD [13] among many others [4, 14]. Previous 
neoclassical transport calculations of the ambipolar electric field at TJ-II include ECH 
plasmas [15, 16, 17] and also medium-density NBI plasmas [18]. In Refs. [16, 18], 
the calculations were compared with Heavy Ion Beam Probe (HIBP) and passive 
emission spectroscopy measurements: qualitative agreement was obtained. A number 
of additional HIBP measurements exist for ECH plasmas [19, 20, 21], and more 
recently for NBI plasmas [22]. Near the edge, the electric field has been measured 
by means of reflectometry [23, 11] and, very recently, by studying mode rotation 
velocities [24]. 
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Ref. [15] includes a comprehensive study of the transport coefficients and the flux 
balances for two ECH plasmas (although no multiple roots were found, see below), 
but an analogous work is missing for medium-density NBI plasmas, where only the 
ambipolar electric field has been shown. Furthermore, lithium wall coating of TJ-II 
has recently allowed [10] transitions to regimes of relatively high density in NBI plas- 
mas, and these plasmas have not yet been described from the neoclassical transport 
point of view. Additionally, many of the effects reported above show a dependence on 
the magnetic configuration [25, 11]. No qualitative changes are expected in the neo- 
classical radial transport of these configurations since the main Fourier components 
of the magnetic field strength do not change too much [26, 15]. Still, the flexibility 
of TJ-II allows for exploring a large set of configurations (on a shot-to-shot basis [27] 
or continuously [28]) and a general study of the variation of the transport coefficients 
and the neoclassical balance may be of interest. More precise calculations for selected 
discharges arc underway 

In this work, we aim to complete the previous neoclassical transport studies in 
view of the recent upgrades in TJ-II operation: we study a low-density ECH plasma 
and discuss, for the first time in detail from the neoclassical transport point of view, 
the issue of multiple roots at TJ-II. We also study a high-density NBI plasma and 
show for the first time the radial fluxes. Finally, we extend the calculations in these 
two plasmas to seven other configurations usually operated at TJ-II. The paper is 
organized as follows: the basic theory is reviewed in Section 2. The monoenergetic 
radial transport coefficient for the 100_44_64 magnetic configuration, the most usually 
operated at TJ-II, is briefly described in Section 3.1. Then, convolution and solution 
of the ambipolar equation yield the radial fluxes and the radial electric field for the 
two plasmas. A Monte-Carlo technique for error propagation allows us to account for 
the convergence problems of DKES [29] for the long- mean- free-path (Imfp), enhanced 
by the complexity of the magnetic configuration of TJ-II. This is done in Section 3.2. 
Finally we have explored part of the set of magnetic configurations of TJ-II with 
DKES calculations. In Section 3.3 we compare the main Fourier coefficients describ- 
ing the equilibria, show how the monoenergetic radial transport coefficient depends on 
the configuration and finally discuss the implications on the radial particle balance. 
Section 4 is devoted to the conclusions. 



2. Determination of the neoclassical radial fluxes 

The radial electric field E r may be obtained from the ambipolar condition of the 
neoclassical radial particle fluxes, which in a pure plasma composed of electrons and 
ions reads: 

T e {E I )=T i (E I ). (1) 

The neoclassical fluxes are linear combinations of the density and temperature 
gradients and the radial electric field [1, 2]. For each species b, the flux-surface- 
averaged radial particle flux Tb and radial energy flux Qb are: 
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Here n is the density, and 2\, and Z^e are the temperature and the charge. We consider 
current free operation, £j| = 0. The thermal coefficients L^ k at each radial position 
can be calculated by convolution with a Maxwellian distribution of the monoenergetic 
radial diffusion coefficient: 

Lh(r,n,T u T e ,Er) = -= / dx 2 e^ x 1+2 ^ +s ^ D n ( ri v\ fi) . (4) 
V 71 " Jo 

The integration variable, x = v/v^ h , is the particle velocity normalized by the ther- 
mal velocity. The monoenergetic radial transport coefficient Z?n, which we discuss in 
Section 3.1, depends on the collisionality v* ~vR/vl and the electric field parameter 
f2 = £ , r /(?;i?o). Here, v is the collision frequency, l is the rotational transform, R the 
major radius and Bq is the (0,0) Fourier harmonic of the magnetic field strength in 
Boozer coordinates. 

From Eqs. (1), (2), (3) and (4), it is clear that the determination of the radial 
electric field in helical plasmas is a non-linear problem due to the Lj^(E r ) dependence. 
This may lead to several solutions or roots of the ambipolar condition [32] when solving 
the ambipolar equation by means of root-finding algorithms. When this happens, 
we select among ion root E\ and electron root E° according to a thermodynamic 
condition: minimization of the generalized heat production rate due to neoclassical 
transport, see e.g. [17, 33]. 7<0 (/ >0) means that the electron (ion) root is selected, 
where 

E° 

i=l r (r i -r c )d J B r . (5) 

Note that the transition region between roots has been imposed to have zero radial 
width. Another option would have been to solve a diffusion equation for E T , instead 
of Eq. (1). 

For error estimate, we follow the Monte-Carlo method described in Rcf. [34]: we 
start with a database of the radial diffusion coefficient Dn(r,v* with the corre- 
sponding error bars. For every value of r, v* and f2, we give a numerical value to flu 
in Eq. (4) by employing a Gaussian random number and then solve Eq. (1). By re- 
peating this procedure a number of times, we obtain averages and standard deviations 
of E T and the other relevant quantities. 

It must be noted that the neoclassical ordering may be violated at TJ-II [35, 36] 
if the widths of the ion drift-orbits are large. This makes the diffusive picture fail, 
and a convective term corresponding to ripple-trapped particles should be added to 
Y\ in Eq. (1), leading to a smaller electric field. The higher the collisionality, the lower 
the correction. The incomprcssibility of the E x B drift, one of the approximations in 
the calculation [37], is valid for the NBI plasmas, where the electric field is small in 
absolute value, but it may fail when the electron root is realized and the electric field is 
close to resonance values. This may also lead to underestimate the ion flux, and thus 
to overestimate the positive radial electric field. Finally, although the monoenergetic 
calculations do not conserve momentum, momentum-correction is negligible for the 
radial transport of non-quasisymmctric stcllarators [39, 38]. 
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3. Calculation and results 

3.1. The magnetic configuration and the radial diffusion monoenergetic coefficient 

TJ-II is a four-period (N = 4) flexible heliac of medium size, with strong helical 
variation of the magnetic axis and magnetic surfaces with bean-shaped cross-section 
and small Shafranov shift. The so-called 100_44_64 configuration (the most often 
employed during TJ-II operation) has a major radius of R= 1.504m, its minor radius 
is a — 0.192 m and its volume-averaged magnetic field is 0.957 T. We stick to the 
vacuum equilibrium calculated using VMEC [40], which has a fiat iota profile, with 
t(0) = -1.551 and t(o) = -1.650. 

The radial diffusion coefficient has been calculated with DKES, in an independent 
simulation for each radial position and value of v* and f2. The distribution function 
has been described with up to 150 Legendre polynomials and 2548 Fourier modes. 
For the description of each magnetic surface, the largest 50 Fourier modes have been 
employed. 

Fig. 1 shows the normalized radial diffusion coefficient D\\=D\\jD\ x calculated 
for several values of the collisionality and the normalized electric field, at p = 0.7. 
The normalization D X1 is the value of Du in the plateau regime for the equivalent 
axisymmctric tokamak, as in Rcf. [41]. 

The Pfirsch-Schhiter (PS), plateau and long- mean- free-path (Imfp) regimes [2] are 
clearly visible in Fig. 1, with a qualitative dependence on v* and Q equal to that of the 
classical stellarator. The same behaviour has been reported for similar configurations 
of TJ-II and other stellarators [15, 41]. This text-book dependence of D xl with v* and 
f2 can be summarized as follows: 

• v regime for high collisionalities, with reduction of transport for large radial 
electric fields due to E x B poloidal precession [42]. 

• Independence from v* and f2 for intermediate collisionality. 

• 1/v dependence for low collisionalities and O = 0, suppressed by the electric field, 
which leads to v and yjv regimes [3]. 

The large error bars for collisionalities lower than 10~ 4 might lead to inaccurate 
results for low-density plasmas, due to the \ jv dependence of the coefficient. In this 
work, these error bars have been propagated to the final results following the method 
described in Ref. [34]: as we will see in Section 3.2, the collisionality is usually high 
enough so that this procedure yields accurate results. 

A database of r>\ x has been built in the (p,^*,0)-space, with p between 0.1 and 
1, v* between 3 x 10~ 5 and 3 x 10 2 and Q between and 1 x 10 _1 . The convolution 
of Eq. (4) requires interpolation and extrapolation in this three-dimensional database. 
The interpolation is done by means of 3-point Lagrange, with u* , £1 and D* x 
in logarithmic scale. Since is calculated at several tens of radial positions, 
interpolation in p is not necessary (but note that Eq. (4) is local, so that radial 
interpolation would not be required at this step). We have made use of the 
asymptotic collisional and collisionlcss limits [41] where possible. Neither the choice 
of interpolation algorithm nor the extrapolation procedure affect significantly the final 
results [43] . Integration has been made by means of Gauss-Laguerre of order 64 (order 
200 yields compatible results). 
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Figure 1. Monoenergetic radial diffusion coefficient at p = 0.7 as a function of 
the collisionality for several values of the normalized electric field. 
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Figure 2. Plasma profiles for the low-density plasma: density and temperatures 
(left); collisionalitios (right). 



3.2. Radial balances 

We calculate the ambipolar electric field for two regimes of TJ-II: a low-density ECH 
plasma (Figs. 2, 3, 4, 5 and 6) a high-density NBI plasma (Figs. 7, 8, 9 and 10). Since 
the walls are coated with lithium [44], TJ-II plasmas have low impurity content, and 
the effective charge may be taken to be Z e g = l everywhere. 

First of all, the two plasmas show some common features: near the edge, both 
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Figure 3. Thermal transport coefficients for the low-density plasma. 




Figure 4. Radial balance for the low-density plasma: ambipolar radial electric 
field (left); ambipolar particle and energy fluxes (right). 



the temperatures and densities drop to zero in a way that the collisionality rises 
some orders of magnitude. As a result of this, the thermal coefficients will be very 
small and the neoclassical transport will be negligible: in this region, transport 
will be completely anomalous. Therefore, the calculated E v may be relevant only 
if the dominant turbulence is electrostatic, and hence automatically ambipolar; the 
neoclassical radial fluxes are shown for the sake of completeness. 

For electrons, one has almost everywhere L\ 2 L\ x ^> L\ 1 , so the temperature 
gradient acts usually as the main drive for both the radial particle and energy flux. 
For ions, the three coefficients Lj k are somewhat closer, and the ion temperature is 
rather flat, so also the density gradient and the radial electric field are responsible for 
the radial fluxes. 

The low-density plasma has a hollow density profile and a peaked electron 
temperature profile, with a central value around 1 keV, see Fig. 2. In these conditions, 
the electrons are in the Imfp regime, except near the edge, where the temperature 
drops. The ion temperature is lower, about 100 cV, and therefore the ions are in the 
plateau regime. 

The thermal transport coefficients are larger for the electrons for p < 0.6, see 
Fig. 3, and so is the temperature gradient. This yields a positive ambipolar radial 
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Figure 5. Radial particle fluxes as a function of the radial electric field at 4 radial 
positions of the ECH plasma: p = 0.22 (top left), p = 0.61 (top right), p = 0.65 
(bottom left) and p = 0.74 (bottom right). 



electric field in Fig. 4. Fig. 5 shows the solution of the ambipolar equation at p = 0.22 
(top left) and p = 0.61 (top right). At both positions, the radial electron flux, which 
is larger, is reduced by a non-zero electric field via the dependence L^ k (E r ). The ion 
radial flux is driven by the convective term E x jT\ (although at ,0 = 0.61, the effect of a 
radial electric field reducing the orbit width is clearly visible). At p = 0.22, the poloidal 
resonance may be playing a role in Y x for E T ~ 1 kV/m: the diffusion coefficient has a 
peak for the value of E r such that the poloidal E x B drift and the poloidal parallel 
velocity cancel out, and then decreases. At outer positions, the electron collisionality 
is higher and the density and ion temperature gradients rise. In these conditions, 
the ion particle flux is larger than the electron flux for p > 0.61. Around p = 0.65 
(bottom left of Fig. 5), two stable solutions for the ambipolar condition appear: as 
in inner positions, a positive electric field (electron root) which drives convectively 
the ions. But also a negative electric field (ion root) is able to bring ion transport 
to the electron level, via reduction of the radial excursions of trapped ions. At these 
conditions, electrons are the rate- controlling species (i.e., the ion radial flux is reduced 
to the electron level), and the particle flux is driven by the density and electron 
temperature gradients. Finally, closer to the edge, the electron transport drops, and 
only a negative electric field is able to restore ambipolarity. Bottom right of Fig. 5 
shows this happening at p = 0.74. The radial dependence of the ambipolar equation 
has been summarized in Fig. 6, where we show the radial particle fluxes as a function 
of p and E x . The intersections of T[(p, E T ) and T e (p,E r ) give the profile of the radial 
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Figure 6. Radial profile of the solution of the ambipolarity equation for the ECH 
plasma. 

electric field E I (p). 

This general behaviour is consistent with that obtained in calculations for 
similar magnetic configurations of TJ-II [15, 17] and to the experimental data from 
HIBP [16, 22]. 

For a wide range of radial positions, 0.3 < p < 0.7, the radial fluxes are very 
slowly decreasing (almost constant) functions of p. This happens because the ions, 
which are the rate controlling species, are driven by a radial electric field with low 
shear, E T ~ dT e /dr. A jump to higher fluxes may be observed for p> 0.7 in Fig. 4, 
and thus a minimum of the radial fluxes at p = 0.7, the point of maximum density 
and temperature gradients. The jump can also be seen in the thermal coefficients, see 
Fig. 3: it corresponds to the change of electron to ion root and it reflects the difference 
in the drift-orbit size for the two different absolute values of the radial electric field. 
The jump is thus larger for the ions, which is consistent with the data shown in Fig. 5 
(bottom left). If E r were obtained by solving a diffusion equation, the transition from 
electron to ion root would be smoother. Yet, the results of the calculations in Rcf. [17] 
suggest that the minima might remain, both in the particle flux and in the energy 
flux; the formation a particle transport barrier at p = 0.7 has been observed in low- 
density transitions in ECH plasmas of TJ-II (see Ref. [21] and references therein) but 
no simultaneous accumulation of energy has been detected. 

During this low-density transition, a double poloidal- velocity shear layer has been 
measured near the edge [11, 24]. Although the formation of this layer is caused by 
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anomalous transport, once the profiles are set the static electric field may be discussed 
in terms of neoclassical fluxes: the electron temperature (and thus the thermal trans- 
port coefficients) drops much faster than the ion temperature in the range 0.7 < p < 0.9, 
and this requires a negative shear to maintain ambipolarity; the opposite happens for 
p>0.9, hence the positive shear in our calculation. 

For NBI plasmas, the electron temperature is flatter and lower, since the plasma 
has reached the ECH-cutoff density. The ion temperature is still flat, and slightly 
higher due to the NBI heating. The lithium wall coating has allowed [10] transitions 
from plasmas of medium density ((n) ~ 2 x 10 19 m~ 3 ) to plasmas of relatively high 
density ((n) ~ 5x 10 19 m -3 ). The studied plasmas are shown in Fig. 7: both the ions 
and the electrons are in the plateau regime, the ions being slightly more collisional. In 
these conditions, the thermal transport coefficients shown in Fig. 8 are larger for the 
ions, whose drift-orbit size is much larger. Since the electron temperature is rather flat 
except near p = 0.7, the plasma is in the ion root everywhere, see Fig. 9. The solution 
of the ambipolar equation for these plasmas is shown in Fig. 10. The situation is 
similar to that of the ECH plasma near the edge: the radial ion flux, being larger, is 
reduced via L\ 2 (E V ) to the electron level, which is in turn determined by the density 
and electron temperature gradients. The radial fluxes are maximum where the density 
and ion temperature gradients are larger. 

The ambipolar E T for the is qualitatively similar to that obtained for a 
similar magnetic configurations (and lower density) in Ref. [18] and to HIBP 
measurements [22]. 

At TJ-II, the ion temperature is usually measured by a charge-exchange neutral 
particle analyzer [46] and the profiles, obtained on a shot-to-shot basis, are not always 
compatible with the ones deduced from spectroscopy measurements [18]. It is therefore 
meaningful to allow for variations in the ion temperature profile in the neoclassical 
transport calculations. If T; were slightly higher and more peaked than the T\ employed 
so far, no qualitative effects would be expected in ECH plasmas, since E r is determined 
by the electron temperature. The ambipolar flux, being the ions the rate-controlling 
species, would probably be slightly larger. For the NBI plasmas, a more negative 
electric field is expected, maybe slightly peaked near the core region, with no major 
changes in the ambipolar flux. 





Figure 9. Radial balance for the high-density plasma: ambipolar radial electric 
field (left); ambipolar particle and energy fluxes (right). 



Configuration name 


(B(T)) vol 


V(m 3 ) 


t (0) 


l{o) 


100_32_60 


1.087 


0.934 


-1.423 


-1.517 


100_38_62 


0.971 


1.031 


-1.492 


-1.593 


100_40_63 


0.960 


1.043 


-1.510 


-1.609 


100_42_63 


0.931 


1.079 


-1.534 


-1.630 


100_44_64 


0.962 


1.098 


-1.551 


-1.650 


100_46_65 


0.903 


1.092 


-1.575 


-1.676 


100_50_65 


0.962 


1.082 


-1.614 


-1.704 


100_55_67 


0.964 


1.073 


-1.659 


-1.739 



Table 1 . Main parameters of the configuration scan. 
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Figure 10. Radial profile of the solution of the ambipolarity equation for the 
NBI plasma. 

3.3. Configuration dependence of the radial fluxes 

Finally, we make a scan in a relevant part the configuration space of TJ-II. We study 
eight magnetic configurations usually explored in regular operation of TJ-II. Their 
main global parameters are shown in Table 1. The rotational transform at the mag- 
netic axis varies about a 20% during the scan (and so does at the edge, since the profile 
shape is kept unchanged) . Since the volume is approximately constant as well and the 
main Fourier terms in the description of the magnetic field do not change too much, 
no large differences in the neoclassical radial transport are expected [26]. Indeed the 
qualitative behaviour is identical to that of Fig. 1. Wc thus make a discussion on 
the radial profile of several quantities that parametrize the dependence on v* of the 
monoenergetic coefficient for the 100_44_64 configuration and the others. 

The two main contributions to the radial diffusion should come from the helical 
and toroidal (&i,o) curvatures and the toroidal mirror (&o,i) m the Fourier 
decomposition of B given by: 

oo oo 

B(p,4>,6)/B {p)= J2 ^bnMcosfaB-Nrufr), (6) 

n— — oo m— 

where <f> and 8 are Boozer coordinates. They are shown in Fig. 11 for three of the con- 
figurations of the scan: 100_44_64, 100_32_60, 100_55_67 (the center and the extremes 
of the scan). 
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Figure 11. Main Fourier components of the 100_44_64 configuration (top left); 
helical curvature (top right), toroidal curvature (bottom left) and toroidal mirror 
(bottom right) for the configurations 100.32.60, 100.44.64, 100.55.67. 



In the PS regime, for fl = , we have oc k x v* [42]. Here, k oc (6i,o) 2 is 
called toroidal curvature. The small correction for large Q depends weakly on t and 
R. 

According to Fig. 11, we can expect the collisional transport of these 
configurations to be quite similar. Nevertheless, modulation of the toroidal mirror 
term &o.i allows for optimization of radial transport in elongated configurations (see 
e.g. Ref [41]). In Ref. [47], an effective toroidal curvature was defined, including bi : o 
and 6o,i, in order to account for this effect. The latter Fourier term has a large relative 
(although not absolute) variation in this configuration scan [34]. From the slope of 
D\i for high v* in Fig. 1 one can calculate an effective curvature K e g defined by: 

D^*^oo,Q = 0) = ^— . (7) 

Therefore, from the slope of for high v* in Fig. 1 one can extract the local value 
of K e ff- The profile of this quantity is shown in Fig. 12. There is a t -2 scaling (see 
Table 1) at every radial position in our set of configurations, so we focus, from now 
on, on the extremes of the scan. Since D\ x oc (ti?g) _1 and v* oc the contribution 
of the PS regime to the thermal transport coefficient is Ln oc (Kefft 2 ^ 2 )" 1 - These re- 
sults altogether yield a neoclassical transport in the configuration 100.32.60 reduced 
with respect to the others (100.44.64 and 100.55.65) via the trivial scaling LhocBq 2 . 




Figure 12. Profile of the effective curvature for all the configurations of the 
study. 



The results are different in the Imfp regime. In the limit case of classical stcllarator 
(i.e., if only 6 lj0 and 61,-1 are non-zero in Eq. (6)), we have £> xl oc (fr^-i) 3 / 2 /^* f° r 
r2 = [3]. If more Fourier terms are non-zero, particle-trapping in local minima of the 
magnetic field leads to enhanced radial transport. However, one can still describe the 
radial diffusion in the Imfp regime in terms of an effective helical ripple e e g [48, 49, 41]. 
This quantity is defined by: 

DlM^0,n = Q)=(±) 2 ^^, (8) 



^3tt, 

obtained from data such as Fig. 1. It contains information of the helical ripple 61,-1, 

as well as of all the other terms in Eq. (6). In Fig. 13 we show the radial profile of the 

effective ripple for three of the configurations: 100_32_60, 100_44_64, 100_55.67. The 

configuration 100_32_60 has a smaller effective ripple, while those of the 100_44_64 

and 100_55_67 configurations are quite close. Since, for the Imfp regime, one has 
3/2 9 

Ln oce c ff -Bq , the latter configurations will have similar radial neoclassical transport. 
Configuration 100_32_60 will have considerably smaller transport, and the reduction 
will be larger than the L\\(xBq 2 of the PS regime. 

These results altogether lead to a reduced radial neoclassical transport for the 
100_32_60 configuration with respect to the 100_44_64, 100_55_67 configurations. The 
reduction is larger for the electrons (which have a larger contribution of the Imfp 
regime) than for the ions (which are more collisional) . Consequently, this will affect, 




Figure 13. Profile of the effective ripple for the configurations 100_32_60, 
100.44.64, 100.55.67. 



via the ambipolar condition, the neoclassical radial electric field, which will be smaller. 
This is consistent with impurity poloidal rotation measurements for ECH plasmas in 
Ref. [18]. The results are shown in Fig. 14 for the profiles of the ECH plasma, see 
Fig. 2: the differences in E r are larger for p<0.7, where the collisionality is lower. The 
radial particle fluxes are mainly reduced for p>0.7, where the ion-root is realized and 
thus the electrons are the rate-controlling species. The results for the NBI plasmas, 
shown in Fig. 15, are similar to that of the ECH plasma for p>0.7: a small reduction 
of the electric field (which now becomes more negative) together with a reduction of 
the ambipolar particle flux. 

4. Conclusions 

We have presented calculations of the radial diffusion monocncrgctic transport 
coefficient for several configurations of TJ-II. By convolution of these coefficients, 
we have calculated the ambipolar radial electric field and the fluxes for two different 
plasmas corresponding to the 100.44.64 configuration. The convolution included data 
with large error bars due to the poor convergence of DKES in the Imfp regime, but 
a Monte Carlo method for error propagation has allowed us to show that the results 
are accurate even for the ECH plasma. 

The results for the ECH plasma are similar to the ones obtained in previous 
works. No such calculations exist for high-density NBI plasmas, but the some of the 
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results shown here stay in qualitative agreement with the experiment. Future work 
includes comparison with HIBP and CXRS measurements for this kind of discharges. 

Small quantitative and no qualitative differences have been found between 
configurations, since the Fourier spectra are very similar. The results predict that 
configurations with reduced L lead, for the same plasma profiles, to slightly smaller 
radial electric field and to slightly improved particle confinement. 

The results shown here extend the knowledge of neoclassical transport and radial 
electric field to regimes and configurations not explored previously at TJ-II. 
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